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The tunneling probability for a system modelling macroscopic quantum tunneling is computed. 
We consider an open quantum system with one degree of freedom consisting of a particle trapped 
in a cubic potential interacting with an environment characterized by a dissipative and a diffusion 
parameter. A representation based on the energy eigenfunctions of the isolated system, i. e. the 
system uncoupled to the environment, is used to write the dynamical master equation for the re- 
duced Wigner function of the open quantum system. This equation becomes very simple in that 
, representation. The use of the WKB approximation for the eigenfunctions suggests an approxima- 

tion which allows an analytic computation of the tunneling rate, in this real time formalism, when 
the system is initially trapped in the false ground state. We find that the decoherence produced 
by the environment suppresses tunneling in agreement with results in other macroscopic quantum 
systems with different potentials. We confront our analytical predictions with an experiment where 
the escape rate from the zero voltage state was measured for a current-biased Josephson junction 
shunted with a resistor. 
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. I. INTRODUCTION 

r-> 

Macroscopic quantum tunneling is a topic of interest that pertains to the boundaries between quantum and classical 
physics. This field has undergone extensive research in recent years as experimental advances have made possible 
the observation of quantum tunneling effects in some macroscopic quantum variables such as the flux quantum 
i-£h ' transitions in a superconducting quantum interference device, or the decay of a zero-voltage state in a current-biased 
Josephson junction 0, 0, 01 0, 0j ■ Macroscopic quantum systems may be modelled by open quantum systems which 
are characterized by a distinguished subsystem, described by suitable degrees of freedom which are subject to physical 
experimentation, within a larger closed quantum system. The degrees of freedom of the remaining system are not 
subject to experimental observation and act as an environment or bath for the distinguished subsystem, which is 
usually referred to as the "system" for short. The environment acts as a source of dissipation and noise for the system 
and produces quantum decoherence. Caldeira and Leggett in two influential papers Q, |(| considered the effect of 
. \ dissipation on the tunneling rate and noted that dissipation always suppresses tunneling; see also Ref. . This work 
was then extended to many other open quantum systems with different system-environment couplings, and different 
potentials for the field [1, S EO, EH H3, EH E3, El ; see Refs. [HG3 for comprehensive reviews. 



Most work on macroscopic quantum tunneling is based on imaginary time formalisms such as the Euclidean func- 
tional techniques which have been introduced in the classical field of noise-activated escape from a metaestable state 
pa, or the instanton approach introduced for quantum mechanical tunneling or for vacuum decay in field theory 
19, 20, 2ll I22L I23L |24|. These techniques are specially suited for equilibrium or near equilibrium situations, but can- 
not be generalized to truly non equilibrium situations. To be able to deal with out of equilibrium situations we need 
a real time formalism which describes the evolution of the quantum system by means of true dynamical equations. 

There are theoretical and practical reasons for a formalism of noncquilibrium macroscopic quantum tunneling. On 
the theoretical side dissipation, for instance, is only truly understood in a dynamical real time formalism. In the 
classical context thermal activation from metaestable states is well understood since Kramers j2^| in terms of the 
dynamical Fokker-Planck transport equation, where the roles of dissipation and noise and their inter-relations are 
known. On the other hand, an open quantum system may be described by a dynamical equation for the reduced 
density matrix, the so-called master equation, or the equivalent equation for the reduced Wigner function which has 
many similarities to the Fokker-Planck equation. However, at present no compelling derivation of the tunneling rate 
is available in this dynamical framework, that might be compared to the instanton approach for equilibrium systems. 
Consequently, the effect of dissipation, noise and decoherence on tunneling and their inter-connections is not yet 
fully understood. On the practical side out of equilibrium macroscopic quantum tunneling is becoming necessary 
to understand arrays of Josephson junctions, or time-de pen dent traps for cold atoms which are proposed for storing 
quantum information in future quantum computers |26L I27L l28t ' 29]. or to understand first order phase transitions in 
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cosmology |3(A l31|. 

In this paper we propose a formulation of macroscopic quantum tunneling based on the reduced Wigner function. 
In recent years we have considered different scenarios in which metaestable quantum systems are described by the 
master equation for the reduced Wigner function. By using techniques similar to those used for thermal activation 
processes on metastable states |25l l32| we were able to compute the contribution to the quantum decay probabi lity 
due to the environment. This was used in some semiclassical cosmological scenarios for noise induced inflation |3J| 
due to the back reaction of the inflaton field, in the context of stochastic semiclassical gravity [j3 13^, |3(| |3?], |38| ; see 
Refs. [39|,|40( for reviews on this subject. It was also used for bubble nucleation in quantum field theory, where the 
system was described by the homogeneous mode of the field of bubble size and the environment was played by the 
inhomogeneous modes of the field |4lLl42j |. and on some simple open quantum systems coupled linearly to a continuum 
of harmonic oscillators at zero temperature |43j . But in all these problems only the contribution to tunneling due to 
activation was considered. The reason is that the non harmonic terms of the potential that induce the pure quantum 
tunneling of the system lead to third and higher order momentum derivatives of the reduced Wigner function in the 
master equation. One has to resort to numerical methods such as those based on matrix continued fractions in order 
to compute decay rates from master equations in this case mm 113. 

In this paper we are able to deal simultaneously with both contributions, namely pure quantum tunneling and 
thermal activation, to the vacuum decay process using the master equation for the reduced Wigner function. The 
computational innovation that makes this possible is the introduction of a representation of the reduced Wigner 
function based on the energy wave functions of the isolated system, i. e. the system not coupled to the environment. 
This representation is useful in a way somewhat analogous to the way the energy representation is useful in the 
Schrodinger equation. It is quite remarkable that in this representation the master equation can be solved analytically 
under certain approximations. The key to this result is that quantum tunneling is already encoded in the energy wave 
functions, which we can compute in a WKB approximation. 

In order to have a working model in a form as simple as possible, but that captures the main physics of the problem, 
we put by hand the effect of the environment with some phenomenological terms that describe noise and dissipation in 
a simple form. It turns out that these terms can be deduced from microscopic physics, when the environment is made 
by an Ohmic distribution of harmonic oscillators weakly coupled in thermal equilibrium at high temperature. Thus 
the model is not strictly suited to describe tunne ling from the false vacuum, or zero temperature transitions. The 
appropriate equations are known in this case |43l |48| and include time dependent noise and dissipation coefficients, 
and anomalous diffusion. Thus the model studied here is a toy model at low temperature. We may expect qualitative 
agreement with relevant experiments at very low temperature, but no precision comparable to the instanton approach 
in this case. We illustrate this by comparing the analytical predictions of our model with an experiment where the 
escape rate from the zero voltage state was measured for a current-biased Josephson junction shunted with a resistor 
0. Since our main purpose here is the introduction of a working formalism for out of equilibrium macroscopic 
quantum tunneling we leave the more realistic and involved computation for future work. 

Master equations play also an important role in elucidating the emergence of classicality in open quantum systems as 
a result of their interaction with an environment. In fact, as the master equation gives the quantum evolution of initial 
states, defined by the reduced Wigner function at some initial time, it has been of great help to study decoherence. 
In particular, it has been used to clarify the decoherence time scales, or the way in which the environment selects a 
small set of states of the system which are relatively stable by this interaction, the so-called pointer states, whereas 
the coherent superposition of the remaining states are rapidly destroyed by decoherence |49t l5Cl 1-5 ll \5'A l53l | . Recently 
the effect of the interaction with the environment on coherent tunneling has been analyzed in the framework of an 
open quantum system that is classically chaotic: a harmonically driven quartic double well |54ll55|. It turns out that 
in this problem, which requires a numerical analysis, tunneling is suppressed as a consequence of the interaction with 
the environment. The dissipation and diffusion terms in the master equation have been derived assuming a high- 
temperature limit of an Ohmic environment, as in the model discussed in the present work. Thus, our calculation can 
also be understood in the context of the study of decoherence on quantum tunneling. In fact, we are able to identify 
the terms responsible for decoherence and their effect on tunneling in our model. 

This paper is organized as follows. In the next two sections ITT1 and ITTT1 we review the theory of tunneling in closed 
systems and introduce the energy representation for Wigner functions. This extended review is necessary both to 
establish our conventions and to recall specific results which are central to the main argument. In section llVl we 
introduce the environment and compute the dynamical equation for the reduced Wigner function of the system. In 
section we introduce the so-called quantum Kramers equation, which is a local approximation to the transport 
equation which captures the essential physics. For the purpose of computing tunneling from the false vacuum state, 
the quantum Kramers equation can be simplified further under the so-called phase shift approximation. In section 
El we a Pply the foregoing to the actual estimate of the tunneling amplitude for the open system. In section IVIII we 
compare the analytical predictions of our model with an experiment on a current-biased Josephson junction. Finally, 
in Section IVIIII we briefly summarize our results. In the Appendixes we provide additional technical details. 
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FIG. 1: A schematic plot of the potential. For an energy E below the barrier there will be three classical turning points, also 
shown. 



II. TUNNELING IN QUANTUM MECHANICS 



In this section we review the WKB method to tunneling in quantum mechanics. The energy eigenfunctions in the 
WKB approximation we obtain will play an important role in the energy representation of the Wigner function that 
will be introduced latter. 



A. The system 



We begin with the simple closed quantum mechanical system formed by a particle of mass M in one dimension 
described by a Hamiltonian 



with a potential U given by 



H-^ + UV), (1) 



U(x) = ±Mn 2 x 2 -^x\ (2) 



for small values of the coordinate x. This is a fairly general potential for a tunneling system, it is the basic element 
in the dashboard potential, which is a very good model for a flux trapped in a superconducting quantum interference 
device (SQUID), or a single Josephson junction biased by a fixed external current pi I5L flfil IHfij . For technical reasons, 
it is convenient to assume that for large x the potential flattens out and takes the value U (x) — —Urx, both negative 
and constant. The tunneling process ought to be independent of the form of the potential this far away from the 
potential barrier. We present a sketch of this potential in Fig. 1. 

There is one classically stable point at x = 0, and one unstable point x — x s = 2MQq/X, corresponding to an 
energy e s = 2M 3 f^/(3A 2 ). The curvature of the potential is U"(0) = MOq at x = 0, and -U"(0) at x s . The other 
point at which U (x) = is x = x ex n = (3/2)x s . For x 3> x ex a the potential flattens out and is constant. 



B. The WKB approximation 



If we assume that the particle is trapped in the potential well, that is in its false ground state or false vacuum, 
the tunneling probability can be computed in this simple problem in many ways. One of the most efficient is the 
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instanton method which reduces to the computation of the "bounce solution" . The most attractive aspect of this 
computation is that it can be easily extended to field theory where the tunneling probability is then interpreted as 
the probability per unit time and volume to nucleate a bubble of the true vacuum phase. The rate for quantum 
tunneling is T c i osec i = (f2o/27r)a g cxp(— Sb/K), where Sb is the action for the "bounce" (or instanton), namely the 
solution to the classical equations of motion which interpolates between x = and x = x ex u in imaginary time 
Sb = 2 Jo"""* dx \J2MU (x), and the prefactor a q = (YIQitSb /K) 1 / 2 . Our expression for the potential is so simple 
that the above integral can be computed explicitly: Ss/ti = 18e s /(5eo)i where Eq = \1t£Iq is the zero point energy of 
a harmonic oscillator with frequency Qq. 

Here, however, we will concentrate on a real time approach by expanding the false vacuum state as a linear 
combination of true eigenstates of the Hamiltonian. To the required accuracy, it is enough to work with the WKB 
approximations to the true eigenfunctions; see for instance Refs. |57L l58j . The instanton method reviewed in the 
previous paragraph can, in fact, be easily justified by this semiclassical approximation. Here we explain in some 
detail this standard procedure to obtain the eigenfunctions by matching the WKB solutions in the different regions 
of the potential. These solutions will play a crucial role in the energy representation for the Wigner functions to be 
introduced latter. 

Let < E < e s be the energy of the particle in the potential well, and ipE the corresponding eigenfunction. The 
Schrodinger equation is 

-ti 2 d 2 

ipE + U (x)iJE = Eip E - (3) 



2M dx 2 
Let us define 

p(x) = y/2M\U(x)-E\, (4) 
and the integral S(x,y) (note the order in the integration limits) 

S(x,y)= f dx'p(x'). (5) 
Jy 

The WKB solutions are obtained from these elements. We have to match the WKB solutions in the different 
regions across the potential function. The details of this calculation are given in Appendix A. The WKB solution 
i/je for energies in the range < E < e s is given by Eq. I|A14|) . where xl < xr < x out are the three classical 
turning points for the cubic potential J2J); see Fig. 1. The normalization constant Ke in Eq. (|A14|I is obtained by 
imposing the continuous normalization of the eigenfunctions given in Eq. l|A"T6)l and it is given in Eq. (|A23(I . Of 
particular relevance is the value of the eigenfunction -0b (x) at values x 3> x ou t- This gives the main contribution to 
the continuous normalization integral. The value of the eigenfunction at i > x ou t, as computed in Appendix A, is 



tPe{x)~J ^ M sin fP^L + S E ) , (6) 
y nnpoo V H ) 

where the phase 8e is introduced in Eqs. (|A24I) and p OQ (E) is defined by Eq. (0J when x ^S> x ou t; see also Eq. <|A17|1 . 

We are interested in the details of the eigenfunctions near the false vacuum state, since we will be dealing with 
tunneling from vacuum. Thus, in the remaining of this section we give explicitly the values of the normalization con- 
stant Ke and the phase shifts Se near this vacuum state. Therefore let us impose the Bohr-Sommerfeld quantization 
condition (|A15|) and let Eq be the corresponding lowest energy, that is, n — in Eq. (|A15|) . As we will see in the 
next subsection this defines the false vacuum energy. Expanding the integral in Eq. (J3J) around Eq we find that close 
to the lowest energy value 

S(x r ,x l )~^-t(E-E ), (7) 

where r is defined by 




(8) 



U (x) - E 



Thus cos 2 (S (xr,xl) /K) ~ (t 2 /H 2 ) (E — Eq) 2 , and evaluating the right hand side of (|A26|) at Eq, we conclude that 
K E has poles at the complex energies 

E± = E Q ± ie, e = A e -2So(x out ,x R )/fi j (9) 
4r 
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which is in agreement with the standard result 58]. To simplify the notation let us call Sq = So (x out , x r) and 
f = f (E ) + it/A, then we have from Eqs. 1A21I) and l|A22[) that the functions A(E) and B(E) for E near E 
are: A (E) = (r/h) exp(S /h)[F^(E) + F + (E)} and B (E) = (-ir/h) exp(S Q /H)[F^(E) - F+(E)], where F-(E) = 
exp(i/o/fi) (E — E-) and F + {E) = exp(— ifo/H) (E — E + ). Notice that neither A nor B vanish at E±. Finally from 
Eq. (|A23|) we can write the normalization constant near the false vacuum energy, as 

2 _ M_ e _ 4Me 2 e 2S °/ h 

E ~^(E-E ) 2 +e*~ {E-EJ){E-E + Y [ ) 

and from Eqs. 1A24|) the phase shifts are 



A5b _ 



2 e 2S /h 



(E - E ) 2 



>E- 


£_ 




E+ 



(A + iB)=e i ^ n J——^. (11) 



Equations ©, (|1(J|> and (|ll|l are the main results of this section. We notice, in particular, the poles of the norm 
and the phase shifts at E± near the false vacuum energy. The strong dependence on the energy of these functions 
near the false ground energy will play an important role in the next sections. 

1. The false vacuum 

Before we start with the computation of the tunneling rate we have to define what we mean by the decaying state, 
all the wave functions we considered so far are true stationary states and, obviously, show no decay whatsoever. We 
need to confine initially the particle into the potential well in its lowest energy. To this end, we introduce an auxiliary 
potential U aux which agrees with U up to x s (where the true potential reaches its maximum value) and increases 
thereafter. We may assume that the growth of l/ aux is as fast as necessary to justify the approximations below; the 
tunneling rate is insensitive to the details of ?7 aux beyond x s . Thus, we define the decaying state ipo as the ground 
state of a particle confined by J7 aux 

It is obvious from the form of the WKB solutions that ipo agrees with ipE U P to x s , i. e. ipo(x) — t[>E (x) for 
x < x s , where Eq is the Bohr-Sommerfeld ground state energy for the auxiliary potential C/ aux , which corresponds to 
n = in the condition l|A15|) . Beyond x s , ipo will decay rapidly to zero, unlike ipE - Like any other wave function, 
admits a development in the complete base of energy eigenfunctions ipE, as 



% (x) = J dE C e ^e (x) , (12) 
where due to our normalization the Fourier coefficients are given by 

C E — J dx tpE {x) ipo (x) . (13) 
To find these coefficients, we observe that ipo (x) is a solution to the Schrodinger equation with the auxiliary potential 



-ft 2 d 2 



2M dx 2 



ipo + U^(x) ip = EqiPo. (14) 



Let us add to both sides of this equation the term [U (x) — ?7 aux (x)]ipo and then multiply both sides by ipE (x) and 
integrate to obtain 

POO 

{E - E ) C E = - / dxip E {x)[U auK (x)-U(x)]ip {x). (15) 

J x s 

An important consideration is that ipo (x) is a smooth function (as opposed to a distribution), and, unlike ipE it is 
normalizable, so Ce must also be smooth. This means that it is allowable to assume E ^ Eq in Eq. (|15|l : Ce can 
then be found by analytical continuation. To estimate the right hand side of Eq. 1)15(1. let us introduce; cf. Eq. 



(x) = V2M|£4ux (x)-Eo\. (16) 
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To the right of x s we may use the WKB approximation with the decaying solution into the forbidden region to write 

ik(x)=il> {x s )e-M:.*«~M dv . (17) 

On the other hand, ipE (x) is given by Eq. (|A12|) in Appendix A. If E is close to Eq, then Eq. J7J applies, and we 
may write 



i\> E (x) ~ 2K t 



T -{E- E ) F + (x s , xr) piy)dy + |f_ (x s ,x R ) K. p[v)dy 



(18) 



Substituting the two previous expressions into the right hand side of Eq. (|15fl we see that we have to compute the 
two following integrals, 

p oo 

J±= dx [C/ aux (x) - U (as)] e - * S:.\p^(v)±P(v)]<hi m (19) 

J x s 

The integral, J_, is dominated by the region near the lower limit, where p mx (x) is close to p (x) and we can write 



from where we obtain 



j_ = h] jf^l-(E-E )j dxe -if: 3 \p^(y)-*>(v)ldy, (20) 

where the remaining integral is made negligible by an appropriate choice of t/ aux . For the other integral, J+, we see 
that the corresponding exponential factor in Eq. (|19|l decays faster than the exponential factor of J_, so that the 
region which effectively contributes to the integral is narrower. Since the pre-exponential factor vanishes at the lower 
limit, we find J + ~ 0. Finally, putting all these pieces together into the right hand side of Eq. (|T5|) we get to leading 
order, 



(E - E ) 

whose solution, assumed smooth, is 



211 (r ) 

C E + 2K E ^o (x s ) t\I -j^ 1 F + (x s ,x R ) 



= 0, 



Ce = -2Ke^o(xs)t^/^^-F + (x s ,x r ). (21) 

We note that Ce is independent of the choice of U aux beyond x s , as it should. 

Thus, we have found the false vacuum wave function in terms of the energy wave functions of the original problem. 
The false ground state is a superposition of energy eigenstates which are fine tuned in such a way as to produce 
destructive interference outside the potential well. Notice that Ce, because of the factor Ke in Eq. I|21|l . peaks near 
the energy of the false ground state, and has a strong dependence on the energy near this ground state energy. 

2. Tunneling from the false vacuum 

Let us now compute the tunneling rate assuming that the particle is described initially by the false ground state 
■00- At time t, we have 



(x, t)= J dE er lEt l h C E i>E (x) , (22) 

The persistence amplitude is 

p(t) = J dx Vo (x) ip Or, t) = J dE e~ lEt/h C 2 E . (23) 
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To perform the integration we can close the contour of integration in the complex E plane adding an arc at infinity, 
whereby we pick up the pole E_ in K\.\ cf. Eq. H21[). Therefore p(t) goes like 



p(t) ~ exp 



4r 



1 exp (~S {x out ,x R/ 



(24) 



provided t is not too large. The tunneling rate for this closed system, r c ; oset j, may be defined from the persistence 
probability p 2 (t) ~ exp(— T c i osec it), so that T c i ose d — (1/2t) exp(— 2So/h), which agrees with the result of the bounce 
solution. Note that if we take the classical lowest energy E = 0, then xr = xl =0, x ou t — x ex it, and Sb = 2S(x ex u, 0), 
but 5*o here is the action corresponding to a particle with false vacuum energy Eq, which differs from zero, consequently 
it differs from Sb/2. This difference is accounted for by the prefactor a q in the instanton result. 

An equivalent way of deriving this result is to estimate the integral by a stationary phase approximation. We can 
write the integral of Eq. lj2l?|) as 



p(t)= J dE e F ^G[E] 



where F [E] = —iEt/K + ln(K E ) with K\ given by (|1(J|> . We consider everything else going into G [E] as relatively 
slowly varying. The stationary phase points are the roots of F 1 = 0. When t — > oo, these roots must approach E±. 
Write e. g. E = E- + er. If |r| < 1, then r = (ih/et)[l — (h/2et) + ...]; this approximation is consistent if et/h > 1. 
At the stationary point 



F[E] = — 



E_t 



In [t] + constant. 



The first term accounts for the exponential decay, as E_ = E — ie where e is given by Eq. (J5J), in agreement with 
the previous result 124fl . On the other hand, F" ~ i 2 , so the Gaussian integral over energies contributes a prefactor 
of order t -1 . 



III. WIGNER FUNCTION AND ENERGY REPRESENTATION 



A very useful description of a quantum system is that given by the Wigner function in phase space, which is defined 
by an integral transform of the density matrix [6(il6ll |. The Wigner function for a system described by a wave function 
tp(x) is 



(25) 



where the sign convention is chosen so that a momentum eigenstate tp Po (x) ~ e %VaX l h j \J2-kYl becomes W Po (x,p) 
(\/2ttH)8 (p — po). Moreover, it satisfies 



/ 



dp W(x,p) = \ip (x)\ 2 , 



dx W (x,p) 



e -ipx/H 

dx — ip (x) 



'2nh 



(26) 



and it is normalized so that J J dx dp W(x,p) = 1. Thus the Wigner function is similar in some ways to a distribution 
function in phase space, it is real but, unlike a true distribution function, it is not positive defined; this is a feature 
connected to the quantum nature of the system it describes. 
The Schrodinger equation for the wave function if;, 



-ti 2 d 2 d 



(27) 



translates into a dynamical equation for the Wigner function, which is easily derived. In fact, by taking the time 
derivative of l|25|l. using the Schrodinger equation (|27|l . and integrating by parts we have 



d_ 

m 



W(x,p) 



^-DK-D-^+DkHl)}- 



f dy Apy/h [ 




d - 


J 2ttH \ 


V M J 


dx . 
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For the cubic potential (J2J) we have U (x — y/2) — U (x + y/2) - , , . ,, , « _ ,.. 

y e ipy/h _ —i]iQ pe vy/fi anc i y3 e ipy/n — ifpg^e %pv ' n , we get the equation for the Wigner function 



: ih 3 d%e ipy ' h 



we 



-MQfixy + (\/2)x 2 y + (A/24)y 3 and, noting that 



dt 



W(x,p) 



d p d 



1 ' dp Mdx + 24 dp 3 



W(x,p), 



(28) 



which may be interpreted as a quantum transport equation. The first two terms on the right hand side are just the 
classical Liouville terms for a distribution function, the term with the three momentum derivatives is responsible for 
the quantum tunneling behavior of the Wigner function in our problem. A theorem by Pawula 44] states that a 
transport equation should have up to second order derivatives at most, or else an infinite Kramers-Moyal expansion, 
for non-negative solutions W(x,p, t) to exist. The above equation for the Wigner function circumvents the implications 
of the theorem since it need not be everywhere-positive. Even if we have an everywhere-positive Gaussian Wigner 
function at the initial time, the evolution generated by an equation such as Eq. (|28l) will not keep it everywhere- 
positive. Thus, here we see the essential role played by the non-positivity of the Wigner function in a genuinely 
quantum aspect such as tunneling. 



A. The energy representation 

Given that a wave function ip can be represented in terms of the energy eigenfunctions i\)e, defined by Eq. J3J), as 

ip(x)= ( dE Ce^e (x) , (29) 



we can introduce a corresponding representation for W(x,p) in terms of a base of functions We 1 e 2 (x,p) m phase 
space defined by 

W El E 2 (x,p) = J ^ e^' h (x - |) r E2 («+!)■ (30) 
Then W (x,p) can be written as 

W (x,p) = J dE x dE 2 C Ei e 2 W Ei e 2 (x,p) , (31) 
where, in this case, we have Ce^e 2 = Ce 1 C E2 . On the other hand from the definition of We 1 e 2 (x,p) we can write 

^-W ElE2 (x,p) W E > lE > 2 (x,p) = J [ip El (x - |) ipE 2 (x + I) Vb; (x - |) ^ (x+ |)} , 
where the p integration has been performed. If we now call z = x — y/2, z' = x + y/2; then dxdy — dzdz', and 



J ^W ElE2 (x,p) W E{E ; (x,p) = ^5 (Ex -E[)S(E 2 - E> 2 ) , 



(32) 



which gives the orthogonality properties of the functions We 1 e 2 - This suggests that any Wigner function may be 
written in this basis as 

W (x,p, t) = J dE x dE 2 C Ei e 2 (t) W Ei e 2 (x, P ) . (33) 

We call this the energy representation of the Wigner function. In this representation, the master equation or the 
quantum transport equation l|28(l is very simple 

Ql C E lE2 (t) = Y {El ~ ° ElE2 {t) ' (34) 

as one can easily verify. One might give an alternative derivation of the tunneling rate from this equation, by taking 
the initial condition for the Wigner function which corresponds to the false vacuum. In fact, in the next section we 
will use the energy representation of the Wigner function to compute this rate in a more complex problem involving 
coupling to an environment. Note that the dynamics of the transport equation in the energy representation is trivial 
and the initial condition is given in terms of the coefficients ()21|l which we have already computed. The task would be 
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more difficult starting from the transport equation in phase space, such as Eq. Q28[l. since the third derivative term 
makes the solution of the equation very complicated. One has to resort to methods such as those based on matrix 
continued fractions in order to compute decay rates from master equations for open quantum systems with third order 
derivative terms 0, liH 14(1 147| . We call the attention to the simila rity of this representation to that based in Floquet 
states that one can use when the Hamiltonian is periodic [62 . l63l l64l l65| . The power of this representation will be 
seen in the following sections when we consider our quantum system coupled to an environment. 



IV. THE OPEN QUANTUM SYSTEM 



So far we have considered a simple closed quantum system. From now on we will consider an open quantum system 
by assuming that our system of interest is coupled to an environment. As emphasized by Caldeira and Leggett Q any 
quantum macroscopic system can be modelled by an open quantum system by adjusting the coupling of the system and 
environment variables and by choosing appropriate potentials. One of the main effects of the environment is to induce 
decoherence into the system which is a basic ingredient into the quantum to classical transition 0, Effi IHH fHTL l58l | . 

The standard way in which the environment is introduced is to assume that the system is weakly coupled to a 
continuum set of harmonic oscillators, with a certain frequency distribution. These oscillators represent degrees of 
freedom to which some suitable variable of the quantum system is coupled. One usually further assumes that the 
environment is in thermal equilibrium and that the whole system-environment is described by the direct product of 
the density matrices of the system and the environment at the initial time. The macroscopic quantum system is then 
described by the reduced density matrix, or equivalently, by the reduced Wigner function of the open quantum system. 
This latter function is defined from the system-environment Wigner function after integration of the environment 
variables. 

In order to have a working model in a form as simple as possible, but that captures the main effect of the environment, 
we will assume that the reduced Wigner function, which we still call W(x, p) , satisfies the following dynamical equation, 



U M ~ M Ox + 2A % W +1 d-p{ P + M ° dp 



W(x,p), (35) 



where 7 which has units of inverse time is the dissipation parameter, and a 2 the diffusion coefficient. The last two 
terms of this equation represent the effect of the environment: the first describes the dissipation produced into the 
system and the second is the diffusion or noise term. An interesting limit, the so-called weak dissipation limit, is 
obtained when 7 — ► 0, so that there is no dissipation, but the diffusion coefficient 70- 2 is kept fixed. We will generally 
refer to equation (|35|) as the quantum Kramers equation, or alternatively, as the quantum transport equation. It is 
worth to notice that this equation reduces to a classical Fokker-Planck transport equation when h = 0: it becomes 
Kramer's equation |25ll32| for a statistical system coupled to a thermal bath and has the right stationary solutions. 

This equation can be derived in the high temperature limit [4Sl IfifiL l67l l68l Ififll FfT |. In fact, assuming the so- 
called Ohmic distribution for the frequencies of the harmonic oscillators one obtains that, in this limit, ct 2 = ksT ', 
where is Boltzmann's constant and T the bath temperature. In the low temperature limit, however, the master 
equation for the reduced Wigner function is more involved |43l Esj . It contains time dependent coefficients and an 
anomalous diffusion term of the type vd x d p W(x,p), where v is the anomalous diffusion coefficient. Nevertheless, a 
good approximation to the ct coefficient is given at zero temperature by a 2 ~ ^MIq. For simplicity we will base our 
analysis in that equation even though we are interested in quantum tunneling from vacuum which means that our 
quantum system is at zero temperature. 

Equation is often used to describe the effect of decoherence produced by the diffusion coefficient to study the 
emergence of classical behavior in quantum systems; this is a topic of recent interest; see Ref. 1531 for a review. Of 
particular relevance to our problem is the study of decoherence in quenched phase transitions |7l| , and the effect of 
decoherence in quantum tunneling in quantum chaotic systems pi. l5q . 

The reduced Wigner function W(x,p) describes the quantum state of the open quantum system, and given a 
dynamical variable A(x,p) associated to the system its expectation value in that quantum state is defined by, 



(A(x,p)) = / dxdp A{x,p) W{x,p) . (36) 
Then one can easily prove from Eq. (!.' j5f) that defining, 

N = J dxdp W(x,p), (E) = J dxdp (j^ + u{x)^W(x,p), (37) 
we have N — and (E) = -j((p 2 /M) - Act 2 ). 
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A. Energy representation of the reduced Wigner function 

Let us now use the base of functions in phase space We x e 2 i, x iP)i introduced in Eq. (|30|l . to represent the reduced 
Wigner function W (x,p, t) as in Eq. (|33|) . The previous N and (E) have very simple expressions in the energy 
representation: 



N 



J dE C EE (t) , (E)= J dEE C EE (t) • (38) 



To check the last equation we note that / dxdp [(p 2 /2M) + U (x)] We 1 e 2 ( x ,p) — E\5 {E\ — E 2 ) , which can be easily 
proved by explicit substitution of the definition of We ± e 2 , an( i trading powers of p by derivatives with respect to y 
into expressions l|30[l . and partial integrations. 

The quantum transport equation l|35(l in the energy representation becomes, 



§- t C El E 2 (t) = ^~{E 1 - E 2 ) C Ei e 2 (t)+jj dE[dE' 2 Qe^e^Ce-^ (t) , (39) 
where, after one integration by parts, 



which has the contributions from the dissipative and the diffusion or noise parts, respectively, as 

Qe 1 E 2 ,E[E' 2 = Q [ E^E 2 .E'^ + Q E^^E'^ ( 41 ) 

From Eq. ll.'iOH it is easy to see that these coefficients can all be written in terms of the following matrix elements: 

X ElE2 = J dxx i\) Ex (x) ^e 2 {x) , (42) 

P Ei e 2 = j J dx ipEi (x) ( x ) i (43) 

( XP ) Ei e 2 = \\ dx x^ Ei (x)^e 2 (x), (44) 

X E X E 2 = J dx X 2 IpEi (x) ips 2 (x) . (45) 

Explicitly, we have that 

QeX.e^ = ^ [( XP )e iE[ & (E 2 - E> 2 ) - P ElE[ X E2E , - X ElE ,P E2E , + (XP) E2E , S (E 1 - E>)] , (46) 



Q 



(N) 



k B T r 



2X Ei e' 1 Xe 2 e! 2 - X% iE ^5{E 2 - E 2 ) - X% 2E ,J(Ei - E[) 



(47) 



Thus, in terms of the coefficients Ce 1 e 2 the dynamics of the quantum transport equation is very simple. This 
equation, in fact, resembles a similar equation when a Floquet basis of states are used |62l l63l 163. 165) . which are very 
useful when the Hamiltonian of the system is periodic in time. The Floquet basis is discrete in such a case and a 
numerical evaluation of the corresponding matrix elements (|42|l - (|45|l can be performed; see for instance |54l |55| for a 
recent application. It is remarkable that in our case approximated analytic expressions for these matrix elements can 
be found. 



B. Some properties of the matrix elements 

The matrix elements (|42|) - (|45|l have a clear physical interpretation and several relations can be derived among 
them. Note that Xe^e 2 is the matrix element of the position operator X in the energy representation. Since 
Xip E (x) = xt[) E (x) , we must have J dE\ X EEl t[) El (x) = xip B (x) . 



11 



On the other hand, Pe 1 e 2 is the matrix element for the momentum operator. The canonical commutation relation 
[P, X] — —iH, implies [H,X] = (—ih/M)P, and taking matrix elements on both sides we have 

{E x - E 2 )X ElE2 = -^P El E 2 - (48) 
Also, X% B2 is the matrix element of X 2 , therefore 

X 2 = f dEX ElE X E E 2 . (49) 



On the other hand, (XP) E E is the matrix element of XP, consequently [(XP) E E ] * = — (XP) E E cor- 
responds to PX, and (XP)^ + {XP) EaEi = [X,P] EiE2 = iM(E x -E 2 ). Also \xP) EiE2 - (XP) E \ Ei = 
(iM/K) (Ei — E-i) X EiE2 , where the commutator [H, X 2 } has been used in the last step, therefore 

( xp ) Ei e 2 = l -^-( E i- E2 ) X kE 2 + j$ (Ei - E 2 ) . (50) 

We have, also, that (XP) EiE2 = J dE X Ei eP E e 2 - One may check, for consistency, that these relations imply N = 0. 
In Appendix B a test of the quantum transport equation in the energy representation (and of the above matrix clement 
properties) is given by checking that a stationary solution with a thermal spectrum is, indeed, a solution in the high 
temperature limit. 



C. Computing the matrix elements 



The matrix elements contain singular parts coming from the integrals over the unbound region beyond x s . These 
singular parts are easy to compute, since far enough the wave functions assume the simple form ©. When performing 
the calculation of the singular parts of the matrix elements we will use that when x — > oo, we have the identities 



sin(px/K) 



np 



Hp), 



cos(px/h) 



0. 



(51) 



which can be easily checked by taking the Fourier transforms of these functions with respect to p. 

The computation of the singular parts of the matrix elements I|42() - (I45|I may be reduced to the evaluation of three 
basic integrals. These integrals are 



Aa,s (Pi,P2) = J dx sin[(pi =pp 2 ) x/% + Si T , 



B{pi,pi)= I dx sm.{pxx/K+ S±) sm(p 2 x/h + <5 2 ) . 



and 



where, for simplicity, we have written pi = p^ (Ei) and Si = S (Ei) (i — 1, 2). The matrix element Xe x e 2 is 

2M 



(52) 



(53) 



X 



E1E2 



Yi-k^PxPi 
M 



dx x sin (pix/h + Si) sin (p2x/h + 5 2 ) 



dA dA 



dpi dp2 \dpi dp 2 



dSi dS 2 



B 



Ky/PlP2 

where A = (A s - A A )/2 and A = (A s + A A )/2. The matrix element X EiE2 is 

dx x 2 sin (pix/h + Si) sin (p2x/h + 82) 



(54) 



X 2 

E\E 2 



2M 



hn^pip 2 

2M 
k^JPiPi 



dC (dSi 



dpi V 9pi 



D 



(55) 
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where, it is easy to show that C = {dB/dpi) — {d5\/dp\)A, and that D = —(dA/dpi) — {dSi/dpi)B. The matrix 
element Pe^e 2 is 

P Bi b 2 ~ - (pi +p 2 A 56 

which according to the relations among matrix elements derived in the previous subsection is related to Xe ± e 2 by 
Eq. I|48[l. The remaining matrix element (XP)e 1 e 2 , on the other hand, can be computed from the element X^ iE2 
according to Eq. ifBll)) 

1. The integrals A(pi,p2) and B{p\,p2) 

Thus, we are finally left with the computation of the integrals (|52|) and (|53|l . The integral B (pi,P2) of Eq. I|53|) is 
dominated by its upper limit x 



1 f s 

B(pi,P2) ~ ^ / cos[(pi -p 2 )x/h + 5i - 62] 



sin [{pi — p 2 ) x/h + Si — 82 



2 (pi - P2) 

- Y*(px-P3), (57) 

The integrals Aa,s {pi,P2) are more subtle. The integral As is clearly regular on the diagonal. Since we are interested 
mostly on the singular behavior of the matrix elements, we can approximate A$ ~ 0. On the other hand Aa is 
exactly zero on the diagonal. Close to the diagonal, the integral is dominated by the region where the argument of 
the trigonometric function is small, and thereby the integrand is non oscillatory. Estimating the upper limit of this 
region as x ~ h {pi — P2)~ , we get 

A A - ft" 1 {pi - pa) x 2 + {Si -S 2 )x= hPV ( — - — J + ... , (58) 

\P1 -P2) 

where the dots stand for regular terms. Actually, this argument would allow us to introduce an undetermined 
coefficient in front of the principal value PV, but in the next section we show that h is the correct coefficient, as 
follows from the canonical commutation relations. 

Thus, we are now in the position to give the explicit expressions for the singular parts of the matrix elements and 
write, finally, the quantum transport equation in its explicit form. This is done in detail in the next section. However, 
there is an approximation we can use that drastically simplifies the computations, and is discussed afterwards, in 
subsection IV CI 

V. THE QUANTUM KRAMERS EQUATION 

In this Section we explicitly compute the quantum transport equation (|35|) satisfied by the reduced Wigner function 
in the energy representation. 

A. Matrix elements 

First, we need to compute the matrix elements described in section llV CI We begin with the matrix element Xe 1 e 2 
which according to l|54|) and H57 M 58 [) can be written as: 



Mh 

^-E 1 E 2 



it dpi \P1-P2J dpi 



\/P\P2 

We go next to the matrix element Pe 1 e 2 ^ which from Q56f) and l|58|l can be written as, 



(59) 



P Ei e 2 = -^L^-{pi+P2)PV—!—. (60) 
2vr pi - p 2 
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These two operators X and P are connected through Eq. I|48|) . It is easy to check that the two previous results satisfy 
this relation. Just notice that from Eq. (|A17|I we can write E\ — E 2 = {p\ — p 2 )/2M which together with Eq. 1|59|) 
for Xe 1 e 2 l ea< i to —ih/M times the right hand side of Eq. H6()[). that is 

(Ei - E 2 )X Ei e 2 = ~j^Pe 1 e 2 - 
Another check of the previous results is the consistency with the canonical commutation relations 

J dE {P Ei eX EE2 - X EiE Pee 2 ) = -ih5 (E l - E 2 ) . (61) 

This check requires a little more work. First it is convenient to change to momentum variables and write, 8 (E\ — E 2 ) = 
(M / y/pip 2 )S (p\ — p 2 ) . Then one needs to compute the integral 



dpPV 



Pi-P 



PV 



P-P2 



-Hir 2 8 (pi - p 2 ) , 



(62) 



The evaluation of this integral is easily performed using the following representation of the principal value 

pv (!p)=j^ ipUh ^^\), 

which is easily proved by taking the Fourier transform of PV(l/p). With the result of Eq. I|62|l it is straightforward to 
check that the commutation relation (|61|l is an identity within our approximation. This consistency check is important 
because it can be used to fix to h the coefficient in front of the principal value of Aa in the argument leading to Eq. 

(EH- 

We can now move to the matrix elements for X 2 . Having an expression for X Ei e 2 in Eq. Ij59(l it is best to compute 
X Ei g 2 directly from the relation l|49|l which leads to 



X 



E1E2 



Mh 2 
VP1P2 
'dth 
,dpi 



1 fdS 1 dS 2 \ d 



op\dp 2 7r \dpi dp 2 J dp 2 



5 {.Pi - P2) + 



Pi -P2 



(63) 



where we have used the result (16211 and performed the E integration or, more precisely, the p integration. 

The matrix element (XP) Ei e 2 = J dEX E\E Pe e 2 can be analogously obtained from the expressions l|59"|) and ||BDJ. 
The result is 



(XP) 



EiE 2 



iMh 
^VpTpi 



d 1 d8 

2^275— $(P1 -P2) + -TT 1 {Pi +P2)PV 

dpi 7r dpi 



1 



Pi - P2 



(64) 



A further consistency check of these expressions comes from the property H5U|) . which is satisfied within our approxi- 
mation. 



B. The quantum transport equation 

Finally, we can write the quantum transport equation in the energy representation l|39(l with the coefficient Q given 
by 141|l . The values of the dissipative and noise parts are given, respectively, by (j46(l and H47JI , which can be directly 
computed using the matrix elements deduced in the previous subsection. It is convenient to define, 

C ElE2 (t) = ~^=C V1V2 it) , (65) 

\JP\P2 

and the result is the rather cumbersome expression (|C1|) given in Appendix |0 As explained there we can get a local 
approximation of the quantum transport equation IjCll) , namely 

dC PlP2 -i , 2 2 , 7 ( d , d 



dt 2Mh iP ^ P2)C ^ + 2{dp-i + dp-2 ] [iPl +P2)C ^ ] 
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It is suggestive to give interpretations to the last three terms in this quantum Kramers equation. The first, of 
course, is the dissipation term, whereas the second and third are diffusion terms. The first involves the dissipation 
coefficient, that defines a time scale tr ~ 7 -1 , which is the relaxation time. 

Before we go on with the interpretation of the different terms, it is important to recall the meaning of the coefficients 
CpiP2 j or C'E 1 e 2 ■ First, we note that these coefficients are directly related to the coefficients Ce of the energy eigenfunc- 
tions which make the tunneling state from the false vacuum in the isolated system, i. e. when there is no interaction to 
the environment. Thus, the coefficients Ce x e 2 clearly give the quantum correlations between wave functions of differ- 
ent energies that make the tunneling system. These coefficients are initially separable Ce 1 e 2 (0) — ^Ei(0)C E (0). In 
the isolated closed system its time evolution, as given by Eq. is simply CE 1 E 2 (t) = Ce 1 e 2 {Q) exp[—i(Ei — E 2 )t/h], 
which means that these correlations keep their amplitude in its dynamical evolution. 

This is very different in the open quantum system. The negative last term in Eq. I|66[) has no effect when E\ = E 2 , 
i. e. for the diagonal coefficients, but its effect is very important for the off-diagonal coefficients. In fact, the amplitude 
of the off-diagonal coefficients exponentially decays in time, on a time scale of the order of 

(67) 

where tr is the relaxation time, \b — h/(2a\fM) is a characteristic de Brolie wavelength (in the high temperature 
case when a 2 = k B T it corresponds to the thermal de Broglie wavelength), and Id ~ a 2 hy/Ea + C^/(eVM) is a 
characteristic length of the problem with a a dimensionless parameter that measures the scale of the energy differences 
of the off-diagonal coefficient, E\ — E 2 ~ as; so it is of order 1 when the energy differences are of order e. The time 
scale ijoTj) can be estimated by taking the derivatives of the phase shifts Si (i = 1,2) near the false vacuum energy 
i?o, which is where the energy wavefunctions pile up. Thus, the last term of equation lj66[) destroys the quantum 
correlations of the energy eigenfunctions. The time scale td may be considered as a decoherence time |49j . and thus 
the effect on tunneling of this term may be associated to the effect of decoherence. 

Another time scale in the problem is, of course, the tuneling time which according to (|24(l and © is given by 
Ttunn ~ hje. Its relation to td is given by Td ~ Ttunn/ (& 4 D), where the dimensionless parameter D is defined in 1|77JI . 

The last of the diffusion terms is the only one that survives in the phase shift approximation which we introduce 
in the next subsection. This is justified by the strong dependence of the phase shifts Si, (i = 1, 2) on the energy near 
the vacuum energy, see Eq. (|11|) , which make the derivatives of these functions very large near Eq . Note also that 
the range of energies (and momenta) in Eq. Ij66(l is also limited to near Eo as the coefficients Ce that describe the 
tunneling state are peaked there; cf. Eq. (|21|l . 




C. The phase shift approximation 



The phase shift approximation is based on the observation made in Section lll Bl that the phase shifts are fast varying 
functions of energy near resonance. This suggests: (a) we only keep the singular terms, and of these, (b) only those 
which contain derivatives of the phase shifts. Under this approximation we can go back to Eqs. H54fl - I|56|) to write, 



Xe x e 2 



-HM fdSi 
Pi \dpi 



5 (pi - P%) 



(68) 



Ei E2 



^P\P2 \OP\! 



P2) 



moreover Pe 1 e 2 ~ 0, and (XP)e 1 e 2 



0. Now we see from gnj) that Q {D) 



(N) 



E\ E2,E-^ E 2 



(69) 



0, so that only the diffusion 



, matters, and this term is proportional to S (Ei — E[) 5 (E 2 — E' 2 ) (we change from variables Pi to 



term Q 

Ei according to pidpi = MdEi). This means that we are working in a weak dissipation limit. Finally, the quantum 
transport equation i|39|) can be written as 



d_ 

dt 



Ce x e 2 (t) — —L [Ei, E 2 ] Ce x e 2 (t) , 



(70) 



where 



L[E 1 ,E 2 ]= l -{E 1 -E 2 ) + ^ f 



Pi 



d8 x 



P2 



dS 2 
dE-2 
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This equation, of course, can also be obtained from Eq. (|66|l in the limit where only the phase shift terms of the 
environment are kept and we return to the E^ variables instead of the pi. In the spirit of the phase shift approximation, 
we shall replace pi (i = 1, 2) by their values at resonance, whereby 

L [E U E 2 ] = i( El -E 2 ) + 2 1 a 2 (E + U*,) (J^- - ||) , (71) 

where, from Eq. Q11JI. the phase shifts derivatives are 

85 _ -i ( 1 1 
~dE ~ ~2 \E-E- ~ E-E- 




VI. TUNNELING IN THE OPEN QUANTUM SYSTEM 



We can now compute the tunneling rate from the false vacuum for our open quantum system. Thus, let us assume 
that our particle at t = is trapped into the well of the potential J5J in the false ground state with the energy Eq, 
i. e. the ground state of the auxiliary potential Z7 aux introduced in Section fllB II We know from that section that the 
wave function ipo of this state can be expressed in terms of the eigenfunctions ipE by Eq. I|12|) with the coefficients 
Ce given by Eq. J3TJ. In terms of the reduced Wigner function, which we may call Wo(x,p), this state is easily 
described in the energy representation l|3*3")l by the coefficients Ce 1 e 2 (0) = Ce 1 {0)C e (0), where Ce(0) is just given 
by Eq. (|21|1 . Because the dynamics of the quantum transport equation is trivial in the energy representation (|70|l the 
time dependence of the coefficients CE 1 E 2 (t) is simply 

C Ei e 2 (t) = exp (-L [Ek, Ei] t) C ElE2 (0) , (73) 

so that, according to Eq. (|33l) . the Wigner function at any time is 

W (x,p, t)= J dE Y dE 2 exp (-L [E U E 2 ] t) C El (0) Ce 2 (0) W Ei e 2 (x,p) . (74) 

From this we can compute, in particular, the probability of finding the particle at the false vacuum at any time. 
In terms of the false vacuum Wigner function and the Wigner function of the tunneling system we may define that 
probability as 

p 2 (t) = 2iTh J dxd P W (x,p)W(x,p,t). (75) 

This equation is justified by observing that in the closed system of Section [H] where the state is described by the 
wave function i/j of Eq. I|22|l and the false vacuum is described by the wave function ip of Eq. I)12H . the square of the 
persistence amplitude l|23|) is given, in fact, by Eq. (|75|l when the definition of the Wigner function , i. e. Eq. (|25|l . is 
used. For the open system the quantum state is not described by a pure state and, in general, the Wigner function 
W(x,p, t) can be written as W = Y^iPiWi where pi is the probability of finding the system in the state fa and Wi is 
the Wigner function for the state 4>i- The definition (|75|) leads in this case to p 2 (t) — '}2 li Pi\{4'o\4 l i)\ 2 ^ which is indeed 
the probability of finding the system in the state ipo. Eq. (|75|l when the energy representation is used becomes 



P 2 



(t) = J dE x dE 2 exp (-L [E U E 2 ] t) C 2 Ei (0) C% 2 (0) . (76) 



To compute p (t) we shall use the stationary phase approximation. The idea is that the integration paths for E\ 
and E 2 may be deformed simultaneously in such a way that the integrand comes to be dominated by Gaussian peaks. 
For late times it is enough to seek the stationary points of L [E\,E 2 ] ■ In principle, we could include K Ei and K E2 as 
fast varying components of the integrand, but these functions are really fast varying in the vicinity of E_ and E + , 
which, when 7^0, are essential singularities of the integrand and must be avoided. Note that when deforming the 
path of integration, we should avoid regions where Re L < 0. Then, calling 

E i = E +er i (z = 1,2), D = jfw 2 {E ° \ Uoo) , (77) 
where e was introduced in Eq. ©, and using Eq. (|72|) for the phase shift derivatives we find that 

te 



L[E u E 2 ]t = 



i(n-r 2 ) + 2D 



l + rj 1 + r, 



(78) 
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The stationary phase condition, dL/dEi = and dLjdE^ = 0, reads 



h [ri,r 2 ] 



8D 



1 



(1 



(79) 



and there is a similar equation / 2 [ri,r 2 ] = 0, where / 2 [ri,r 2 ] is obtained from /i [ri,r 2 ], with the substitution 

of the multiplying factor ri/ (l + rf) by r 2 / (l + r|) . Observe that / 2 [ri,r 2 ] = — (fx [^2 ' r i ] )* ^ so a solution of 
Eq. 179(1 with r 2 = r\ is automatically also a solution of / 2 [ri,r 2 ] = 0. We shall seek stationary points of this 

kind. It is clear that for a solution of this kind, r-yj (l + r 2 ) must be real. So, writing r\ — £ — iry, we must have 



4£ 2 (i + e 2 ) = (1 + e 2 



7? 2 ) , and solving for £ 2 we find 



3e = \(2-v 2 ) 2 



1/2 



-1-rf*. 



(80) 



We notice that for each value of ij there will be two possible solutions for £. The path of integration must go through 
both of them. Substituting, n = £ — irj, into the complex stationary phase equation Q79JI we obtain two equations for 
£ and ij. The real part of Eq. (|79"|) leads to the previous Eq. ljSU|) . which is independent of D, and the imaginary part 
leads to 



V 



4?? 2 (1 



V 



if 



4£V 



For each value of D, we are interested in the solutions with the lowest possible positive value of r\. 
Finally, the contribution of each saddle point to the integral will be 

4(i)=A ± e-*™, 



(81) 



(82) 



where 



2e?7 



1 - WD^tj 



e-v 2 ) 2 



4£V 



and where the prefactor A± depends on the second derivatives of L. Comparing to the persistence probability 
/^closed ~ ex P {—2te/h) for the isolated closed quantum system, which follows from the persistence amplitude Eq. l(2~l*|) 
and the definition of e given in Eq. @, we conclude that the ratio R — T ' open /T 'dosed of the tunneling rates between 
the open, T open = L, and closed, T c i osec i — 2e/H, systems is 



R = r)[l- 16DS_ 2 r] (1 + £ 



V 



(83) 



where the parametres £ and 77 are solutions of the algebraic equations l|8U|) and (|81|) with the lowest possible positive 
value of rj. 

In Fig. 2 we plot R as a function of D from a numerical solution of Eqs. I|81|l and l|8U|l (full line); we see that R < 1 
for all values of D, with tunneling being strongly suppressed when D is large. In the asymptotic limits D — > and 

1 /3 

D — > 00 it is possible to obtain analytical approximations. In the former, 77 = 1 — <5, £ = ±<5, where S = (D/2) . In 
this limit 



Rd- 



2 V 2 



1/3 



(84) 



which is also plotted in Fig. 2 (dashed line). In the opposite limit, D — > 00, we have r) = 8/ {27 D) , £ = ±1/V3 and 

4 



Rd^c 



27 D' 



(85) 



which is plotted in Fig. 2 (dotted line). 

Observe that if we included the factors of K\. {i — 1, 2) into the fast varying part of the integral, then the saddle 
points shift by an amount ~ h/{et); therefore these results are reliable when {et/fijD 1 / 3 3> 1. On the other hand, 
when et/H is extremely large, the energy integrals are dominated by the contribution from the lower limit — and 
the decay rate turns to a power law. 
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FIG. 2: The ratio _R between the tunneling rates for the open system and the closed system as a function of the dimensionless 
parameter D (full line). We also plot the analytic approximations for low and large values of D (dashed and dotted lines, 
respectively) . 



In this section we confront with experimental results. One of the experiments in which macroscopic quantum 
tunneling has been observed in more detail is a single Josephson junction between two superconducting electrodes 
biased by an external current. The macroscopic variable in this case is the phase difference (p of the Cooper pair wave 
function across the junction [1|. A macroscopic system always interacts with an environment and a physical Josephson 
junction is generally described by an ideal one shunted by a resistance and a capacitance which phenomenologically 
account for the effects of the environment. 

The basic equations relate the (gauge invariant) phase difference <p, the current / and the voltage V across an 
ideal junction j.ii l56| - A zero voltage supercurrent I = I c simp, should flow between two superconducting electrodes 
separated by a thin insulating barrier, where I c is the critical current, the maximum supercurrent the junction 
can support. If a voltage difference V is maintained across the junction the phase difference evolves according to 
2eV — hdtp/dt, where 2e is the charge of a Cooper pair, and fi/(2e) is the flux quantum. It is convenient to introduce 
the characteristic energy scale Ej = HI c /2e; note that the work necessary to raise the phase difference across the 
ideal junction from to <p is (— Ej) cos ip. A real junction is modelled as an ideal one in parallel with an ordinary 
resistance Jin, which build in dissipation in the finite voltage regime, and an ordinary capacitance C, which reflects 
the geometric shunting capacitance between the two superconducting electrodes. The so-called bias current flowing 
through the device is I = I c simp + V/Rn + CdV/dt, and substituting the previous relationship between V and dip/dt, 
this equation becomes a differential equation for the phase difference of the Cooper pair tp: 



which has the form of the equation of motion for a particle in a one-dimensional potential with friction. We have 
introduced the "mass" M = h 2 C / (2e) 2 , the friction coefficient 7 = 1/(RqC), and the "potential" 



where s = I/I c . This is the so-called "tilted washboard" model. Note that when s = 1 the local minima of the tilted 
cosine become inflection points, and that no classical stable equilibrium points exist when s > 1. We will be interested 
in the case in which s is smaller but close to 1, that is when the external biased current I is slightly less than the 
critical current I c . 

In this case, the potential may be approximated by a cubic potential in the neighborhood of any stable stationary 
point. Let us consider the stable stationary point closest to tp = 0. The stationarity condition is simp = s and the 



VII. COMPARING WITH EXPERIMENT 




(86) 



u i<f) = (~ e j) ( cos f + sip) , 



(87) 
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stability condition is cosp > 0. We see that there must be a solution ipo < ir/2, thus let us write ipo — f — k, where 
k ~ y2 (1 — s) ~ \/l — s 2 . We henceforth introduce a new variable a; = — <^o and the shifted potential 

U (x) = U ((pa + x) - U (ipo) , (88) 

which leads to the same equation IjSfijl than the potential (|87|l . In the new variable, the stable stationary point lies 
at x = 0. The closest unstable stationary point, the maximum of the potential lies at x s — 2k, and following the 
notation of Sec. [HI we have that the height of the potential barrier has an energy e s = U (x s ), given by 

s s ~\ej [2(1 - s)f 2 ~ 2 - Ej (l-s 2 f 2 . (89) 

Finally, the next root of the potential is x ex u, where U (tpo + x ex u) — U (tpo), which can be approximately written 
as sx ex i t — sm(x ex u — k) ~ k. Observe that x ex u ~ 3k ~ (3/2)x s , which is similar to what happens with the cubic 
potential of Sec. [H] Note that when s is nearly 1 the height of the potential barrier e s is much less than the potential 
difference between adjacent wells and the potential can be approximated by a cubic potential. For |a;| no much 
larger than x ex u, we may approximate U (x) ~ e s (x/x s ) 2 (l — x/x ex i t )/(l — x s /x ex it), from where we may define the 
frequency f^o of small oscillations around x = 0: Mfi 2 , = U"(0), which gives 

2 3e s Ae 2 Ej , 2 ,i/ 2 _ 2 

n °~2(l-^) M ~-^C-( 1 - a ) =U " (90) 

where lo p is the "plasma frequency" of the junction. We may introduce lu p0 = (Ae 2 Ej /f^C) 1 ^ 2 — (2eI c /HC) 1 ^ 2 . 

In the literature there are several reported observations of tunneling in this or more complex set-ups 0, 0, 0, [H 
l73l 173. FH 17^. l77l[78l . From this wealth of data we have chosen to focus on comparison against the experiments 
reported in Ref. j4|. In that paper, tunneling was observed under a fixed bias current, as opposed to a time-dependent 
one |72j |. The fixed current environment is closest to the ideal situation to which our model applies. 

A critical input in comparing theory to observations is the value of the critical current for the junction. In Ref. 
the critical current is not derived from the tunneling rate itself, as advocated in Ref. 80], but extracted from an 
independent set of measurements at high temperature. We therefore find that the comparison of theory and experiment 
may be done in two ways: it is possible to derive the escape temperature from the independently measured critical 
current, as it is done in Ref. or else it is possible to induce the value of the critical current from the observed 
rate. This value may then be compared to the one obtained in high temperature determinations. In the following we 
present both analysis. 



A. Tunneling rates for the open system 

Let us first compare our theoretical model against the value for the tunneling amplitude (extrapolated up to) at 
T = 0, as reported in Ref. 0. The relevant values for the bias current /, the critical current I c , its ratio s, the 
self-capacitance C, and shunt resistance Rq are: 

I - 24.710 fiA, I c - 24.873 fj,A, s - 0.9934, C - 4.28 pF, R Q - 9.3 Q. (91) 

These parameters were measured in the classical limit. In previous similar experiments on current-biased junctions 
but with low dissipation 0> Q it was shown that at low temperature the tunneling rate became almost independent 
of temperature. 

Since h = 1.054572xl0~ 34 Js and e = 1.602176x 10~ 19 C we have, in the conditions of Ref. 0], Ej - 8.185xl0- 21 J, 
which using Boltzmann's constant ks = 1.380650 x 10~ 23 JK _1 , may be converted into a temperature Ej/ks ~ 
592.9 K. Therefore the height of the barrier e s , uj p o, the plasma frequency f2o, the zero point energy £0 = hHCto, an d 
the friction coefficient 7 in our model are: 

e s /k B - 589.74 mK, lu p0 - 132.88 x 10 9 s" 1 , Q ~ 44.918 x 10 9 s - \ 

e /k B - 171.55 mK, 7 - 25.123 x 10 9 s _1 . (92) 

We observe that the zero point energy is lower than the barrier, although not much lower. In fact, there is only one 
trapped excited state. 

It is often convenient to give tunneling rates in terms of an effective escape temperature T esc . This escape temper- 
ature is defined Q from a given tunneling rate, T, by the equation 



(93) 
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For the closed system, either the WKB approximation or the instanton method yield T esc ~ 71 mK, with a barrier 
penetrability A, defined by T c i OS ed = (l/2r) exp(— A) (which gives A = (2/K)Sq in the WKB approximation), of 
A = 8.459; see Eq. l|D7|) in Appendix D. 

For the open system, the tunneling rate in our model may be expressed in terms of i?, introduced in Eq. (|83[) as 
the ratio between the open and closed systems rates, namely 

R 

I" open — -R r closed — ~z ^ • (^^) 

It 

The relevant parameter in R is D which was introduced in (|77|l . where at zero temperature tr 2 ~ eo and also £o ~ Eq. 
Thus we have D — 167r 3 (7/J7o)(l + C^oo/£ ; o) exp(3A). With the numerical values given in l|91() - (|92() . D is very large and 
we are in the limit of Eq. ijSEfr. namely R ~ 4/ (27D). Comparing Eq. with Eq. i|9"3|) . and assuming U^IEq < 1, 
we may find our predicted effective escape temperature from e s / (fcsT esc ) ~ A + ln(27Z?/4) ~ 4A + 8. 827, which leads 
to T esc ~ 14. 255 mK. 

The experimental result [4| when extrapolated to the zero temperature limit is, in terms of the effective escape 
temperature, T esCtexp ~ 45 mK. The experimental result is in good agreement with the instanton prediction of 
Caldeira and Leggett for weak dissipation at zero temperature p| ; see also 0, ^| . Our prediction goes in the 
sense that dissipation also suppresses tunneling but the predicted value for the effective escape temperature is lower 
than the observed result. We recall that our prediction is based on Eq. (|70|l which was derived under an approximation 
that fully neglects any activation that would increase the effective escape temperature. A better approximation based 
on Eq. i|66[l would require a numerical solution, but our starting phenomenological equation (|35|) is too crude to 
expect a quantitative agreement with the experimental result at very low temperatures at the level of accuracy of the 
experiment. 



B. The critical current 

The value of the critical current is a crucial input in the calculation of the tunneling rate |80j . We are primarily 
interested in the tunneling rate under fixed bias current. However, the critical current was not determined under these 
conditions, but extracted from the average of measurements of tunneling rates under a ramped bias current performed 
at several different temperatures. For this reason, it is meaningful to contrast theory and experiment in a different 
way than in the previous subsection, namely, instead of computing the tunneling rate from the given values of the 
critical and bias currents, we may instead compute the critical current from the given bias current and tunneling rate. 
We may then see if the critical current obtained this way is consistent from that obtained from the ramped current 
measurements also reported in Ref. |4|. 

Recalling that the relation between the critical and bias currents is I c = I/s, where we will now assume that s is 
so far unknown but close to unity, we may write Ej — HI/(2es). Then, from <|89|l we have that the height of the 
potential barrier in terms of the parameter s is, e s ~ (?iJ/3e)s _1 (l — s 2 ) 3 / 2 , and from (|90() the frequency of small 
oscillations is flo ~ (2eI/HC) 1 ' 2 s^ 1 ' 2 (l — s 2 ) 1 / 4 . We shall assume that the ground state energy is close to the zero 
point energy Eq ~ eo = MIq/2. The ratio of the ground state energy to the height of the barrier is thus, in terms of 
s, 

Eq p^ 

E S ~ (l- S 2)5/4' 

where p = (9e 3 /2ftC/) 1 / 2 - 1.288 x 1CT 3 . 

As we discuss in Appendix D, due to the form of the potential the action under the barrier, S(x ou t,Xfi), for the 
energy E is equal to the action inside the barrier, S(xr,xl), for an energy E re f — e s — E. Writing as in Eq. I|D4|) 
E re f = 2e s ( (k re f), where ((k) is defined in Eq. I|D5|) . we get (1 - s 2 ) 5/4 s~ 1/2 = p[l - 2C (fc re /)]~ 1/2 , which relates 
the parameter s with k re f. The corresponding penetrability, defined here as A = (2/K)S{x out ,XR), is 

A = E F ^ (96) 
8 1-2C(fc re /)' 1 j 

where F (k re f) is defined in Eq. (|D6|I . The tunneling rate according to Eq. (|94|l is then 

O 2 4 

r = -^2 ! e -4A (Q7) 

open 2^ 7 27 x 16 xtt 3 ■ 1 ' 
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We must compare this with the experimental tunneling rate which may be computed from Eq. (|93|l with the effective 
escape temperature T esc exp ~ 45 mK, and with the barrier hight e s computed with a critical current I c ~ 24.873 /iA. 
Finally, equating the tunneling rate (|97JI with the experimental one, using the relevant numerical values given previ- 
ously, we get 



This equation must be solved simultaneously with the previous equation relating s with k re f. The solution is k re t ~ 
0.1162 and s ~ 0.9968. Thus the critical current is 



Observe that our model neglects any contribution that leads to activation, and therefore underestimates the tunneling 
rate |l4j |. In matching against experiment, this effect is compensated by lowering the predicted critical current. 
Therefore this result must be regarded as a lower bound. With this in mind, the agreement with the value of 
I c ~ 24.873 /iA extracted from the ramped current measurements [4( is satisfactory. 



To conclude, let us briefly summarize our findings. We have computed analytically the effect of decoherence on 
quantum tunneling in a model representing a quantum particle, trapped in a local minimum of a potential, which 
is coupled to an environment characterized by a dissipative and a diffusion parameter. We have used the master 
equation for the reduced Wigner function, or quantum transport equation, that describes this open quantum system. 

Our computational method involves the introduction of an energy representation of the reduced Wigner function 
which is based on the energy cigcnfunctions of the isolated system. The dynamical equation in this representation, Eq. 
(13911 . is an equation for some coefficients that describe the quantum correlations between eigenf unctions of different 
energies. The equation can be explicitly written in a local approximation, the quantum transport equation l|66() . that 
captures the essential physics of the problem. 

In our problem, where the particle is initially trapped in the false vacuum, the transport equation is dominated 
by a term that destroys the quantum correlations of the cigcnfunctions and is, thus, responsible for decoherence. 
The strength of this term is characterized by the dimensionless parameter D, defined in Eq. <|77[1 . which is directly 
proportional to the energy difference between the false and true vacuum. It does not seem surprising that decoherence 
suppresses tunneling, as it destroys the fine tuning among the energy eigenfunctions that makes tunneling possible 
in the isolated closed system. The analytic result for the ratio of the tunneling rates between the open and isolated 
closed systems is given by Eq. (|83|) . This ratio decreases, implying more suppression, with increasing D. 

Our model is based in the phenomenological transport equation l|35|l for the reduced Wigner function, which is a toy 
model at zero temperature. Besides this assumption, in the paper we work at two different levels of approximation. 
We first derive Eq. (|66l) . which still retains the leading effects of noise, dissipation and decoherence. To be able to 
isolate the effects of decoherence, we then proceed to derive the simpler Eq. (|70|l . The actual predictions for tunneling 
rates are obtained from this later equation, which in practice means to pick up the environment terms that give the 
most decoherence, and presumably the most tunneling suppression. This means that we have neglected from the 
remaining terms any activation that would increase the effective escape temperature. Therefore the results from our 
model, derived from Eq. (7UJ|, must be regarded as a lower bound on the actual rates. This is consistent with our 
goal, which is not to provide an alternative to the instanton calculations in equilibrium, but to propose a starting 
point for a real time formulation of noncquilibrium macroscopic quantum tunneling. When seen under this light, it is 
encouraging that the detailed matching against experimental results shows that our model not only captures the main 
effect, since indeed the measured rate is substantially lower than the quantum prediction for the closed system, but 
also yields a suitable estimate of the critical current from the given tunneling rate. One might try to improve on this 
prediction by numerically solving the quantum Kramers equation l|66|l . on which we expect to report in a separate 
publication. The extent of the discrepancy can be seen also as a check on the validity of the phenomenological terms 
introduced into the master equation for the reduced Wigner function at low temperature. 

One should note that our approximations are valid when the initial state is the false vacuum. The results might 
differ, even qualitatively, when more general initial conditions are assumed and the terms that we have neglected in the 
quantum transport equation become relevant. It may not be possible in such a case to solve analytically the quantum 
transport equation. Yet, we should emphasize that this equation in the energy representation is much simpler than 
in the standard phase space representation. 




(98) 



I c ~ 24.789 /j,A. 



(99) 
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APPENDIX A: WKB SOLUTION 

In this Appendix we solve the WKB problem posed in section III Bl The starting point are Eqs. ©, (0} and (J5J 
with the cubic potential of Eq. we have to match the WKB solutions in the different regions across the potential 
function. 

1. Matching from forbidden to allowed regions 

Let Xq be a classical turning point U (xo) — E, and let U' (.To) < 0. Then to the left of xo we have a forbidden 
region, the two corresponding independent WKB solutions of the Schrodinger equation J2J are 

e ±S(x ,x)/h 

F ± (x Q ,x) = -====, (Al) 
whereas to the right of xq the two independent solutions are 

G±(x,x )= (A2) 

y/2 P (x) /h 

and we wish to find the corresponding matching conditions. For x — > Xq , we can Taylor expand the potential around 
xo and write p (x) = k(xq — x) 1 ^ 2 and S(xo,x) — ^k(xq — x) 3 ^ 2 , where we have introduced k — y/2M \ U' {xq)\. 
Similarly for x — > Xq , we have p (x) — n (x — xq) 1 ^ 2 and S (x, xq) — |k (x — xq) 3 ^ 2 . 

If we write x — xq — e lw (xq — x) then iS (x, xq) — S (xo, x) and it would seem that simple analytical continuation 
yields G + (x,x$) — * e _I7r / 4 F + (xq,x). However, this is impossible, recall that if we define the flux J = —i(ip*d x ?p — 
ipdxip*) then the Schrodinger equation implies flux conservation d x J — 0. Now G+ (x, xo) has J = 1 and therefore it 
cannot turn into F + (xq, x) , which is real, and has J = 0. Thus, we try instead 

G+ (x, xo) -» e-™'^F+ (so, a;) + /3F_ {x Q , x) , (A3) 
and imposing flux conservation we obtain (3 — (1/2) exp(i7r/4). We therefore find the matching conditions 

e^-/4 F+ (a . Q) x) + 1 e ±ftr/4 F _ (xQj x) ^ q ± (a;) Xq) t (A4) 



from were we finally obtain, using Eq. (|A2I 



and 



F + (x Q ,x)^ 7 ==== cos[-S(x,x ) + j), (A5) 



F- (x , x) -» 7=== sin I \s (x, x ) + 7 ) ■ (A6) 



y/2p(x)/h 

2. Matching from allowed to forbidden regions 

Now consider the case when U' (xq) > 0. To the left of Xq, we have an allowed region and the solutions are oscillatory 
G± (xq, x), to the right of the turning point we have a forbidden region and the solutions are a linear combination of 
(jAip . By exactly the same procedure of the previous section, after imposing flux conservation across Xo we obtain: 

G T (x , x) -> e ± ^/ 4 F+ (x, x ) + ^^F_ (x, x ) . (A7) 
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Note from these equations that the solution that matches a decreasing exponential is 

1 /1„, . tt\ 1 



(1 IT 
— ; sin — S (xo, x) + — 

^2p(x)/h \h 4 



-F- (x,x ) ■ 



(A8) 



3. WKB solution for < E < e s 

We can now put all this together to write the energy eigenfunctions for our cubic potential @ for energies in the 
range < E < e s . There are three classical turning points in this case xl < xr < x out . To the left of xl we have a 
forbidden zone extending to -co, so we have 

ipE {x) ~ K E F^ (x L ,x); x < x L , (A9) 

where Ke is a normalization constant to be determined latter. To the right of xl we have from (|A6|) 



which after using the definition JSJ) can be rewritten in the region xl < x < xr as 

i> B (x) ~ K E ^ s{x ^ XL)/h - 7r/4) G- (x R , x) + e -<s(^L)/n-^/i) G+ ^ ^ 



(A10) 



(All) 



This expression is in the form suitable for extension to the forbidden region, that is, to the right of xr. Thus, by 
using (|AT|) we have to the right of xr 



ip E (x) ~ 2K E 
which can be rewritten again as 

i> E (x) ~ 2K E 



l " s ( -^S(x R , x L )j F + (x,x R ) + ~ sin (^rS (x R ,x L ) } F l.c.c/,.) 



(A12) 



cos 



(jS{x R ,x L )j e s ^ x ^' h F_ (x out ,x) + isin (~S(x R ,x L fj e - s ^' x ^ h F + (x out ,x) 



which is in a form suitable for extension to the right of x ou t- 

2K E 



(A13) 



i>E (x) 



2 COS (j_S(xR,X L ^j e S(*ou t ,*n)/H sin (^_ S (x,X out ) + J 



+ ~ sin Qs (x R , x L )) e -s<«..*.**)/» cos Qfi (x, x out ) + J 



Note that if we impose the Bohr-Sommerfeld quantization rule 



S {x R , x L ) = -(1 + 2n)H, 



(A14) 



(A15) 



only the subdominant, exponential decreasing part survives. This would correspond to the case when the far right 
region is forbidden and may be used to define energies for false states trapped into the potential well, in particular 
n = will correspond to the false ground state. 



4. Normalization 



All that remains now is the determination of the normalization constant Ke which can be done from the normal- 
ization of the wave functions. The eigenfunctions are subject to continuous normalization 



dx ipEt (x) i/)e 2 (x) =6(E 1 — E 2 ) ■ 



(A16) 
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Since the functions themselves are regular, the singular behavior must come from the upper limit, see for instance 
[57]]. For large enough x, we have from Eq. (JSJ 



■ Poc = y/2M (E + Uoo). 



Let us write from Eq. I jAlj l. 



S(x,x out )=p 00 {x-x out )+ / dx' y/2M (E~U (x')) - y/2M (E + U^) 



if this integral converges, we may take the upper limit of integration to 00 , whereby 

S (x, x out ) = PooX + f (E) , 
where f(E) stands for the second term of (|A18ll . Then, for x ^> x ou t, we can write from (|A14|) and l|A19(l 



1p E (x) 

where A(E) and B(E) are given by 



2hK f 



A (£) sin (^) +£?(£) cos 



A(E) = 2 cos -S(x R ,x L ))e 



S(x out ,XFt)/h 



cos 



■ l -»m(jS(x H ..xi.))e : " 



(A17) 
(A18) 
(A19) 
(A20) 



(A21) 



B i E ) = 2 cos ( ifi si) J e S(*~*.«*)//> sin (/ (£) + \ 



-isin Qs^xl) ] f - s ' 1 -' '"'cos 



Substituting Eq. ljA"20|) into iJXTBJl, the singular terms in the normalization integral are 

dpi 



(A22) 



dx $ El (x) ^e 2 (as) ~ ?i 2 7r— ^ [^l 2 + B 2 (Si)] 



dEi 



6 {E 1 - E 2 ) . 



where the delta function comes from the x integration which brings S(p\ — p-i) and where we have defined pi = Poo (Ei) 
(i = 1,2) and changed from momentum to energy variables according to Pidpi — AldEf, see Eq. (|A17J| . The 
normalization condition reduces to 1 the coefficient of the delta function above 



h\^[A 2 (E 1 )+B 2 (E 1 )] =1. 



This suggests the introduction of the phase Se as follows, 



, M M 
K E A(E) = x I — C0S 5e, K e B{E) = J— s inS E - 



(A23) 



(A24) 



Thus, the eigenfunction at x ^> x out is Eq. ©, that is 



2M . / Poo x \ 
Sm {— + SE ) 



j hirpoo V h 

To work out the constant Ke in greater detail we note that form Eqs. l|A21jl and (|A22|I we have 

A 2 + B 2 = 4 cos 2 (x R , xl)^ e 2S ( x "^t,x R )/h + 1 gin 2 ^ ^ e -2S(* „ tl x R )A 



(A25) 



which is non vanishing as long as E is real. However, if we allow for complex energies, as is typical of unstable states, 
it may be zero provided 



cos- I -S (Xr, X L ) 



^^(\s{x R ,x L ))< -'<••' " 



(A26) 



The left hand side is zero whenever the energy satisfies the Bohr-Sommerfeld condition (|A15|) . 
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APPENDIX B: THERMAL SPECTRUM 

In this Appendix we check that the quantum transport equation l|39|) admits a stationary solution with a thermal 
spectrum. This can be seen as a test on the restrictions satisfied by the matrix elements 142|) - (|45|l with a 2 = fc^T. 
An unnormalized thermal density matrix in the position representation reads, 



p (x, x') = J dE e-^^E (x) ^ E {x') , (Bl) 

where (3 — (ksT) -1 and its associated Wigner function is 

W (x,p) = J dE e~ pE W EE (x,p), (B2) 

which in the energy representation in the base We 1 e 2 of Eq. (|33|) corresponds to the coefficients Ce 1 e 2 = 
e~ l3El 5 (Ex — £'2). Inserting this into the transport equation we get 

J dE er l3E Q El E 2 ^EE = 0, (B3) 

which after using Eqs. (|41|l . (|46|) and 147|) can be written in operator language as 

= ~(XPe-f >H -e-f ,H PX-Pe-f ,H X+Xe-P H P) 

+ _L (X 2 e-? H + e-P H X 2 - 2Xe-P H X) . 

At the infinite temperature limit, (3 = 0, this is 

= — - (XP - PX) — \ (X 2 H + HX 2 - 2XHX) . (B4) 
M n n 2 ' 

The first term is the commutator which gives — M _1 , and the second term can be written as [X, [X,H]), which 
using [H, X] = (h/i)(P/M) is easily seen to cancel the first term. 



APPENDIX C: QUANTUM TRANSPORT EQUATION 



Here we write explicitly the quantum transport equation l|39f) in the energy representation. The coefficient Q in 
Eq. (|3~9"jl is given by (fTT)! . and the values of the dissipative and noise parts of this coefficient are given, respectively, 
by Eqs. fUjl) and (J37J. These parts can be directly written using the matrix elements deduced in Section [V] When 
the coefficients C PlP2 defined in Eq. lt)5|) are introduced the transport equation becomes, 
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where we have used the shorthand notation P(x) = PV(l/x). 

This equation simplifies considerably if we assume that C piP2 ~ C (pi) S (pi — P2), with C (p) slowly varying. This 
is justified by noticing the effect of the second local a 2 term which is negative. This term has no effect on the 
diagonal terms, when p\ = p2, but its effect on the off-diagonal coefficients is very important. In fact, it exponentially 
reduces the coefficients C PlP2 on a time scale of the decoherence time, as discussed in Section[V] One may argue that 
since tunneling is a long time process with a typical scale of time Ttunn ~ h/e the local approximation should give a 
reasonable approximation to the transport equation l|66|l whenever td <C Tt unn . 

Now, from Eq. (|62[l and the diagonality assumption for the C PlP2 coefficients we have that 



/ 



,2 



This term together with the two local terms involving second order derivatives of the momenta in Eq. IjCljl become 

7 m. 2 ( A + A V c pil 



1P2 ' 



k dpi dp 2 , 

from where Eq. I|66|> follows as the local approximation of the quantum transport equation (|C1|) . 

APPENDIX D: TUNNELING RATES FOR THE CLOSED SYSTEM 

In this Appendix we review the calculation of the quantum mechanical tunneling rate for the closed system, that 
is, ignoring the interaction with the environment. The quantum tunneling rate as given by the instanton calculation 

mum is 



r (inst) _ a q Ao , , 

1 closed — 2 T ' [ ' 



where r = 7r/fi , A = S B /h = 18e s /(5e ) ~ 12.376, and the prefactor a q = (1207rA ) 1/2 - 68.306 . With these 
values, the escape temperature defined in Eq. H93|) is 

Tt st) = ,, T' k , B ^ 72.345 mK. (D2) 

3.6 - {e /e s ) lna g 

It is interesting to check that this result agrees with the result we obtain when the dissipation is zero. We can use 
our WKB result as obtained in Sec. |nl see Eq. I|24|l . to write 

V {WKB) _ J_ a (m) 

1 closed ~ 2t ' 

where A = (2/h)S (x out , xr), with 5*0 defined in Eqs. 0)-©, where the potential U(x) is given by Eq. (JHSJ. 

For a cubic potential, the relationship among the energy E, the frequency Q and the action S (xr, Xl) is best given 
in parametric form, 

E = 2e s ((k), n = Oof(k), S(x R ,x L ) = ^-F(k), (D4) 

"0 

with < k < 1, and 

C(fc) 1 ( 2 , 3 ( 1 + fc2 ) ^ + fc2 ) 3 ] (D5) 
f(k) = l^-[4Q(k)f 4 K[k 2 ] 



5/4 



{a (fc) E [k 2 ] - (1 - k 2 ) b (k) K [k 2 ] } , (D6) 



where E [fc 2 ] and K [fc 2 ] are the complete elliptic integrals, and we have introduced the functions Q (k) = 
(1/4) (1 + 14A: 2 + k 4 ) , a (k) = (16/15) (2 - k 2 f - (1/5) (l - k 2 ) (21 - 5k 2 ) and b (Jfe) = (8/15) (2 - k 2 ) - (l - k 2 ) . 
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The Bohr-Sommerfeld condition Eq. (|A15|) for the ground state (n = 0), corresponds to the parameter kcs such 
that F (has) — tt£o/ £ s which implies that kcs ~ 0.1152. This corresponds to ( (has) ~ 0.1423 and / (has) ~ 0.9550, 
while the harmonic approximation for the potential yields 0.1454 and 1, respectively. 

To compute the barrier penetrability, A = (2/h)So(x ou t,XR), we observe that S (x ou t, xr) at energy E is equal to 
S(xr,xl) at energy E re f — e s — E. The exchange of E by E re f is equivalent to the exchange of k by k re f, where 
C (k ref ) = 1/2 - ((k) . For k GS we obtain k ref ~ 0.2433 and F (k re f) ~ 2.4073. Therefore 

A = — F (fc re/ ) - 8.459. (D7) 

This is to be compared against the instanton exponent Aq — lna 9 ~ 8.152. In terms of the escape temperature, the 
WKB approximation yields 

T^ C KB) = jp~77 , , e °( k * (n JTT-s ~ 70.869 mK, (D8) 
F {kref) ~ (eo/e s ) In (n GS /n ) 

which is in good agreement with the instanton result. This, of course, should not be surprising since for a closed system 
our method reduces to the standard WKB calculation. The purpose of this exercise is just to check the consistency of 
our calculation and to illustrate how the instanton and WKB methods compare. That the difference between exp(Ao) 
and exp(A) is accounted for by the prefactor a q of Eq. i|Dl|) can be seen analytically by a perturbative calculation. 
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